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Abstract 

We present new singular solutions of the biharmonic nonlinear Schrodinger equation 
iipt{t, x) - AV + W^^ip = 0, X G R'', 4/d < cr < 4. 
These solutions collapse with the quasi self-similar ring profile tpQ-^ , where 

IV'QB(t,r)K-J--QBf^— r = |xl, 



4-CT 



L{t) is the ring width that vanishes at singularity, rmax(f) ~ roL°'{t) is the ring radius, and a — cr(d-i) 
The blowup rate of these solutions is for 4/d < a < 4, and slightly faster than 1/4 for ct = 4. These 
solutions are analogous to the ring-type solutions of the nonlinear Schrodinger equation. 

1 Introduction 

The focusing nonlinear Schrodinger equation (NLS) 

zVt(t,x)+AV' + |^P"V' = 0, V(0,x) = Vo(x)ei/i(K'^), (1) 

where x e and A — dxixi + ■•■ + dx^^x^: admits solutions that become singular at a finite time, 
i.e., limi^Tc IIV'IIhi ~ where < t < Tc- Until a few years ago, all known singular solutions of the 
NLS were peak-type. By this, we mean that if assume radial symmetry, and denote the location of maximal 
amplitude by 

- = argmax |'(/;|, r — y + ■ ■ ■ + x'^, 



then r,iiax(i) = for < t < Tc, i.e., the solution peak is attained at r = 0. In recent years, however, new 
singular solutions of the NLS were found, which are ring-type, i.e., rniax(0 > for < t < Tc. 
In this study, we consider the focusing biharmonic nonlinear Schrodinger equation (BNLS) 

zV't(i,x)-A2V'+|V|'"^ = 0, V(0,x) = Vo(x) ei/^M''), (2) 

where A^ is the biharmonic operator. Singular peak-type solutions of the BNLS have been studied in |FIP02[ 
IBFMlOi [BF] . Singular ring- type solutions of the BNLS with ct > 4 were studied in |BFG09| . The goal of 
this work is to find and characterize singular ring- type solutions of the BNLS with 4/(i < ct < 4. 
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1.1 Singular solutions of the nonlinear Schrodinger equation (NLS) - review 

The NLS (HI) is called subcritical if ad < 2. In this case, all solutions exist globally. In contrast, solutions of 
the critical (ad — 2) and supercritical {ad > 2) NLS can become singular at a finite time. 

Until a few years ago, the only known singular NLS solutions were peak-type. In the critical case ad = 2, it 
has been rigorously shown |MR03j that peak- type solutions are self-similar near the singularity, i.e., ip ^ ipfi, 
where 

and r = |x|. The self-similar profile R{p) is the ground state of the standing- wave equation 

R"ip) + ^—^R' -R+ IRI^/'^R = 0. 
P 

Since R attains its global maximum at p = 0, ipR is a peak-type profile. The blowup rate of L(t) is given by 
the loglog law 

In the supercritical case ad > 2, the rigorous theory is far less developed. However, formal calculations 
and numerical simulations |SS99j suggest that peak-type solutions of the supercritical NLS collapse with the 
self-similar V's profile, i.e., ip ~ ips, where 

i^s{t,r) = -^^S{p)e^\ (4a) 
' (4b) 



and S{p) is the zero-Hamiltonian, monotonically-decreasing solution of the nonlinear eigenvalue problem 

S"ip) + ^—^S'-S + i^(-S + pS']+\S\^''S = 0, 5'(0) = 0, (4c) 



p 2 

where k is the eigenvalue. Since |>S'(p)| attains its global maximum at p 0, ips is a peak-type profile. The 
blowup rate of L(t) is a square- root, i.e.. 



L{t) - n^T~t, t ^ Tc, (4d) 

where k > is the eigenvalue of (|4c| . 

In the last few years, new singular solutions of the NLS were discovered, which are ring-type |FGW05| 
IFG W07I [Rap06l [RH09l IFG08] . In particular, in |FGW07| . Fibich, Gavish and Wang showed that the NLS (HD 
with d> 1 and ^ < cr < 2 admits singular ring- type solutions that collapse with the ipq profile, i.e., "0 ~ i^Qi 
where 

i^Q{t,r) = -^Q(p)e-+— (5a) 
'^^ ._'^-^max(t) r^,,(i)^roL"~-(t), (5b) 



L^{-^y ^ m ' 

and 

2-(7 _^ ad~2 
"""^ " a{d-l) " ^ a{d-\y ^^^^ 
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Figure 1: Illustration of ring radius rinax(i) and width L{t). 



The self-similar profile Q attains its global maximum at p = 0. Hence, rinax(^) is the ring radius and L{t) is 
the ring width, see Figure [TJ 

A unique feature of the ipg profile (|5ap is the linear combination of the two radial phase terms. The 
first phase term a^LsIx'''^ describes focusing towards r — 0, and is the manifestation of the shrinking of the 
ring radius r^ax to zero. The second term (1 — q;nls)|x'(^ ~ 'max(i))^ describes focusing towards r = rmax, 
and is the manifestation of the shrinking of the ring width L{t) to zero. The discovery of this "double-lens" 
ansatz was the key stage in the asymptotic analysis of the tpg profile, which enabled the calculation of the 
shrinking rate a^Ls, see ([Sc]) and the blowup rate p, see ([7]). 

The NLS ring-type singular solutions can be classified as follows, see Figure [D 

A. In the subcritical case {ad < 2), all NLS solutions exist globally, hence no singular ring- type solutions 
exist. 

B. The critical case ad = 2 corresponds to Q!nls = 1- Since r^ax{t) ~ roL{t), these solutions undergo an 
equal-rate collapse, i.e., the ring radius goes to zero at the same rate as L{t). The blowup rate of L{t) 
is a square root. 

C. The supercritical case 2/d < a < 2 corresponds to < a^Ls < 1- Therefore, the ring radius r,jiax(i) ~ 
roL"'*'"^(i) decays to zero, but at a slower rate than L{t). The blowup rate of L{t) is 

L{t) ^ k{T, - t)P, (6) 

where 

= 1 = 1 

i + Onls 2 - ^JJ^ 

D. The supercritical case a — 2 corresponds to a^Ls = 0, i.e., lim r^a^U) = ryastxiTc) > 0. Therefore, 

the solution becomes singular on the d-dimensional sphere |x| — rmax(7c), rather than at a point. The 
blowup profile ipq is equal to that of peak- type solutions of the ID critical NLS, and the blowup rate 
is given by the loglog law ([3]). 

E. The case a > 2 also corresponds to a standing ring. The asymptotic profile is not given by tJjq, however, 
but rather by the asymptotic profile of peak-type solutions of the ID supercritical NLS. The blowup 
rate is a square root |BFG09) . 

Thus, NLS ring-type singular solutions are shrinking (i.e., limt-^Ta ''max(^) = 0) for ^ < c < 2 (cases B and 
C), and standing (i.e., < limt_j.j-^ r^-^^xit) < oo) for a >2 (cases D and E). 
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Figure 2: Classification of singular ring- type solutions of the NLS, as a function of a and d. A) subcrit- 
ical case - no singular solutions exist. B) crd = 2: equal-rate ipQ solutions [FGWOS] . C) 2/d < a < 2: 
shrinking i/jq solutions |FGW07I . D) cr = 2: standing iJjq solutions |FGW07| [RapMl IREOQ] . E) ct > 2: 
standing non- ^pQ rings |BFG09j . 



1.2 Singular solutions of the biharmonic NLS - review 

The BNLS equation ([2|) is called subcritical if ad < 4, supercritical if ad > 4, and critical if ad = 4. In the 
critical case, equation ([2|) can be rewritten as 

#t(t,x)-A2^+|^|**/''V = 0, i^(0,x) = Vo(x)ei/2(jRd)^ 
The BNLS conserves the "power" {L^ norm), i.e., 

p(o^p(o), p{t)^m)\\l 

and the Hamiltonian 

H{t) ^ HiO), Him] = IIA^II^ - T^U\\ltX% (9) 
In the radially-symmetric case, the BNLS equation ^ reduces to 

#t(i,r) - A^V' + I^I'^^^O, V(0,r) =7Ao(0, (10) 

where 

is the radial biharmonic operator. In particular, the radially-symmetric critical BNLS is given by 

*7At(t,r) - A^V' + I^I'/'^V^O, 7A(0,r) =7Ao(r-)- (12) 

All solutions of the subcritical BNLS exist globally in [FIP02j . In the critical case, they exist globally 
if the input power is below the critical power: 
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Theorem 1 ( [FIP02] ). Let ||i/'o|l2 < Per' where P^^ — WRgW^, and Rb{p) is the ground state of the standing 
wave eguation 

- A^R^ip) -Rs+ \Rn\^"R, - 0, (13) 
with a ~ A/d. Then, the solution of the critical focusing BNLS ([8]) exists globally. 

Numerical simulations |FIP021 IBFMIO] indicate that solutions of the critical and supercritical BNLS can 
become singular at a finite time, i.e., \mit-^T^ llV'llff^ = oo, where < Tc < oo. At present, however, there is 
no rigorous proof that the BNLS admits singular solutions, whether peak-type or ring-type. 

In [BFMIO] , we rigorously proved that the blowup rate of all singular solutions of the critical BNLS 
is bounded by a quartic root: 

Theorem 2. Let %p he a solution of the critical BNLS ^ that becomes singular at t — T^ < oo, and 
let l{t) = Then, 3K = if (||-0o|l2) > such that 

l{t) < K{Tc - t)i/4, < i < Tc. 

We also proved that all singular solutions are quasi self-similar: 

Theorem 3. Let d > 2, and let tp{t, r) be a solution of the radially- symmetric critical BNLS ilS^) with initial 
conditions ipo{r) G -ff^^adiai' ^^^^ becomes singular at t = Tc < oo. Let l{t) = ||A-!/)||2 and let 

Smt,r)^f'^(tm, Kt)r). 

Then, for any seguence tj. — > Tc, there is a subsequence tk, such that S{'ijj)(tk,r) ^'(r) strongly in L', for 
all q such that 

2<q<oo 2<d<A, 

2<q<^^ 4<d. ^^^> 



Since the L^-norm of S'('(/') is conserved, and the convergence of S{ip) to 4* is in L'^ with q > 2, the solution 
becomes self-similar in the singular region (the collapsing core), but not everywhere. Consequentially, the 
solution has the power-concentration property, whereby a finite amount of power enters the singularity point, 
i.e.. 



where P^^J is the critical power for collapse |BFM10[ ICHK 



Peak-type singular solutions of the critical BNLS (IT2]) were studied asymptotically and numerically 
in |FIP021 IBFMIO) . The asymptotic profile of these solutions is 

^RB{t,r) = -^i?B (^) (15) 



Ld/^{t) ^ \L{t) 

where R^{p) is the ground state of ([T3| . The blowup rate of L{t) is slightly faster than a quartic root, i.e., 

lim^^ = (« (16) 

t^T, (Tc - t)P \ OO p > i ' 



Specifically, in the one-dimensional case, the quasi self-similar profile is 



^Rs^nit,^) = TrirTT^B.iB (t^) e^'^'^ (17a) 
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and i?B,iD is the ground state of 

-R^""{O-Ru + \RufB.^ = 0. (17b) 

Peak- type solutions of the supercritical BNLS (flO)) were studied asymptoticahy and numerically in |BF) . 
The asymptotic profile of these solutions is 

^ i^'*- (lit)) 

where (p) is the zero-Hamiltonian solution of a nonlinear eigenvalue problem 

- S^{p) + (^S^ + psS^ - Ap^B + I^bI'^^b = 0, 5b'(0) - 5b "'(0) = 0, (19) 

and K is the eigenvalue. The blowup rate is exactly p = 1/4, i.e., 

L(t)^K{T,-tf^\ (20) 

where k > is the nonlinear eigenvalue of p^ . 

Ring-type singular solutions of the supercritical BNLS ([T0| with cr > 4 were studied asymptotically and 
numerically in |BFG09| . These solutions are standing rings, i.e., lim rinax(0 > 0. The self-similar profile of 

these standing-ring solutions is 

= ^s«,..(i,x = r - r„,a.(0) = ;^5^^B,iD {^^^^^^^ e^^' (21) 

where i/'Sb id(^7^)j see (|18p . is the profile of the peak- type singular solution of the one-dimensional super- 
critical BNLS with the same value of a. The blowup rate is given by (^0]). 



1.3 Analogy of NLS and BNLS 

Table!??] lists the major findings of the previous works |FIP021 iBFMlOi IBFG09] on singular solutions of the 
BNLS, side by side with their NLS counterparts. In all cases, the results for the BNLS mirror those of the 
NLS exactly, "up to the change of 2 -J> 4" . 

We note that current BNLS theory is still missing a key feature in NLS theory, which is the BNLS 
analogue of the quadratic radial phase terms of the asymptotic profiles. Therefore, our asymptotic analysis 
of the BNLS singular solutions produces weaker results than those of |FGW07| . Hence, in this work we "fill 
in" the missing results by relying on the above analogy of the NLS and BNLS, up to the change 2 — > 4. 



1.4 Summary of results 

In this study, we consider ring- type singular solutions of the BNLS (fTO| with 4/(i < ct < 4. We show 
numerically that such solutions exist, and are of the form V'(^, ?") ^ V'Qb (^i where 

I^QB(i,OI = ;^57^QB(p), (22a) 

P = "—J^^ ^max(i) ^ r,L^{t), (22b) 

and 

4 - cr ad - 4: , 
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NLS 


BNLS 


peak-type solutions; 


critical case [BFMlOi 
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2 


ad = A 


asymptotic profile 




/ r > 

U(t)y 






/ r \ 

\L(t)/ 








blowup rate 


slightly faster than 1/2 


slightly faster 


than 1/4 


power concentration 


yes 


yes 


peak-type solutions: supercritical case iBFl 




ad>2 


ad> 4 


asymptotic profile 










( r \ 
\L{t)) 


i ft 






blowup rate 


^1/2 




= 1/4 




power concentration 


no 


no 


"supercritical" standinaj-ring; solutions |BFG09j 




a>2 


a > A 


asymptotic profile 


same 


as ID peak 


same 


as ID peak 


blowup rate 


same 


as 1_D peak 


same 


as ID peak 



Table 1: A comparison of the properties of singular solutions of the NLS and BNLS. These properties 
are analogous, "up to the change 2 — > 4" . 



The ■(/'Qb profile is the BNLS analogue of the ipQ profile of the NLS. Unlike the ipq profile, however, we do 
not know the expression for the double-lens ansatz of V'Qb ■ 

In Section [21 we consider the case a — 4. In this case a = 0, i.e., the solution is a singular standing 
ring. Informal asymptotic analysis and numerical simulations show that the blowup profile is the self-similar 
profile 

Mt,r;a^4) = ;^i?B.o (^— ^f^) (23) 

where i?B,iD is the ground state of (|12l) with ct = 4 and d — 1, and that the blowup rate is slightly faster 
than a quartic root. In other words, the blowup rate and profile are the same as those of peak-type singular 
solutions of the one-dimensional critical BNLS, see ^TE\\ . 

In Section [3l we consider the case A/d < a < 4, for which < aB < 1, see (|22cp . From power 
conservation we deduce that a > Ub- By analogy with the NLS, we expect that a — a^. Therefore, the ring 
radius rniax(i) ~ rQL°'{t) decays to zero, but at a slower rate than L{t). By analogy with the NLS, we also 
expect that the blowup rate of these ring solutions is given by (O with 

1 1 



Numerical experiments support these predictions. 

In Sectional we consider the critical BNLS {a = A/d), which corresponds to Ob = 1. Since the singular 
part of the solution has to be self-similar in r/L, see Theorem |31 a must be equal to unity. By the analogy 
with the NLS, the blowup rate is conjectured to be 1/4. 

In summary, the BNLS singular ring- type solutions can be classified as follows (see Figure [3]): 

A. In the subcritical case {ad < 4), all BNLS solutions exist globally, hence no collapsing ring solutions 
exist. 

B. The critical case ad = A corresponds to Q!b = 1 (equal-rate collapse). The blowup rate \s p — 1/4. 
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B 

ad = A (critical) 
a^ = I (equal-rate) 



C 

A/d < (7 < 4 
< a. < 1 



A 

ad < A 
subcritical 




1 



Figure 3: Classification of singular ring-type solutions of the BNLS as a function of a and d. A: sub- 
critical case (no singularity). B: critical case, with equal-rate collapse. C: A/d < a < 4, shrinking rings. 
D: cr = 4, standing rings. E: cr > 4, standing rings |BFG09] . 



C. The supercritical case 4/d < cr < 4 corresponds to < Ob < 1, hence the ring radius rniax(i) decays to 
zero, but at a slower rate than L{t). The blowup rate is p = - — Jj_4 . 

D. The case a — 4 corresponds to as — 0. Hence the solution is a singular standing ring. The self-similar 
profile Qb is equal to that of peak-type solutions of the ID critical BNLS, and the blowup rate is 
slightly above p = 1/4. 

E. The case cr > 4 was studied in |BFG09| . In this case, the solutions are of the standing-ring type, 
the self-similar profile is equal to that of peak-type solutions of the ID supercritical BNLS, and the 
blowup-rate is a quartic-root. 

Thus, up to the change 2 — > 4, this classification is, indeed, completely analogous to that of singular 
ring- type solutions of the NLS (see Figure [2]). 

1.5 Numerical Methodology 

The computations of singular BNLS solutions that focus by factors of 10* necessitated the usage of adaptive 
grids. For our simulations we developed a modified version of the Static Grid Redistribution method (RWOOi 
IDG09j . which is much easier to implement in the biharmonic problem, and is easily extended to other 
evolution equations, such as the nonlinear heat and biharmonic nonlinear heat equations jBFGOQ) . The 
method of j DG09) also includes a mechanism for the prevention of under-resolution in the non-singular 
region. We extend this mechanism to prevent under-resolution in the transition layer between the singular 
and non-singular regions. See Section [5] for further details. 
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Figure 4: top: Blowup rate of singular ring solutions of the BNLS. The blowup rate increases mono- 
tonically from p = 1/4 at cr = A/d to p = (1/3)— ai a — 4—. For cr = 4 (full circle) p = 1/4 (with 
a loglog correction?) and for cr > 4, p = 1/4. bottom: The shrinkage parameter of singular ring 
solutions of the BNLS. For A/d < a < A, decreases monotonically from 1 to 0+ (shrinking rings). 
For cr > 4, Ob = (standing rings). 



In Figure m (top) we plot the blowup rate p of singular ring solutions of the BNLS, see ©. As a increases 
from 4/d to 4—, p increases monotonically from ^ to ^ — . At cr = 4, the blowup rate is slightly faster than 
a quartic root, i.e., p ^ \- Finally, p = 1/4 for cr > 4. Since 

lim n = l/3, lim p—-, 

the blowup rate has a discontinuity at cr = 4. 

The above results show that a = A is a critical exponent of singular ring solutions of the BNLS. Intuitively, 
this is because the blowup dynamics changes from a shrinking-ring (cr < 4) to a standing-ring (cr > 4), see 
Figure 2] (bottom). We can understand why cr = 4 is a critical exponent using the following argument. 
Standing-ring solutions are 'equivalent' to singular peak solutions of the one-dimensional NLS with the same 
nonlinearity exponent cr [BFG09| . Since cr = 4 is the critical exponent for singularity formation in the one- 
dimensional NLS, it is also the critical exponent for standing-ring blowup. An analogous picture exists for 
the NLS, wherein the phase transition between standing and shrinking rings occurs at cr = 2 jBFG09) . 
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2 Singular standing rings {a = 4) 



In what follows, we show that collapse of ring-type singular solutions of the BNLS with ct = 4 is "the same" 
as collapse of peak- type singular solutions of the one-dimensional critical BNLS. 

2.1 Informal Analysis 

We consider ring-type singular solutions of the supercritical BNLS ([T0|) with cr = 4, that undergo a quasi 
self-similar collapse with the asymptotic profile 

Mt,r) = -^^QAp)e^'^'' p='—Jj^. (24) 

Here and throughout this paper, by quasi self-similar we mean that '0 ^ "^Qb the singular ring re- 
gion r — Tmax = 0{L) , or p = 0(1), but not for < r < oo. 

The asymptotic profile (j24p describes a standing ring if lim rmaxlO > 0- We expect ring-type singular 

solutions of the BNLS with cr = 4 to collapse as standing rings, for the following two reasons: 

1. By continuity, since ring- type singular solutions of the BNLS with cr > 4 are standing rings |BFG09j . 

2. By analogy with singular ring-type solutions of the NLS with cr = 2, which are standing rings [FGW071 
|Rap06[ [R509l. 

In the ring region r — r,„ax — 0{L), as L — 0, the terms of the radial biharmonic operator (ITT|) behave 

as 



4_fe r 



M k = 4 



Therefore, A'^ip ~ d^tp. Hence, near the singularity, equation ((TO)) reduces to 

V'(i,r)-Vrrrr+|V'lV = 0, 

which is the one-dimensional critical BNLS. Therefore, the singular solutions of the two equations are asymp- 
totically equivalent, i.e., 

C=f,<i(^' ^) ^C=td=S^ x = r - r„iax(i)), 

where ip^^^^^i is a peak-type solution of the one-dimensional critical BNLS. 

The above informal analysis suggests that the blowup dynamics of singular standing-ring solutions of 
the BNLS pH]) with d> 1 and cr = 4 is the same as the blowup dynamics of singular peak solutions of the 
one-dimensional critical BNLS: 

Conjecture 4. Let d> 1 and a = 4, and let ip be a singular ring-type solution of the BNLS (jlOp . Then, 
L The solution is a standing ring, i.e., lim rmax(i) > 0. 

2. In the ring region, the solution approaches the ipQg self-similar profile, see (j24p . 

3. The self-similar profile ijjQg is given by 

i^Qsit, r) = iPrbad {t,X = r - TmaxW) , (25) 

where ^jfi^ ^^{t,x), see (1171) . is the asymptotic profile of the one- dimensional critical BNLS. 

4. Specifically, Qb{p) = RsjoiO' where Rg.ioiO ground state of p7b)) . 

5. The blowup rate of L{t) is slightly faster than a quartic root, see (|16p. 
In Section 12.21 we provide numerical evidence in support of Conjecture |4l 
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Figure 5: A singular standing-ring solution of the supercritical BNLS (fTO|) with d — 2 and cr = 4. A) 
Ring radius r^ax as a function of the focusing level 1/L. B) The rescaled solution, see (^5]) . at L{t) — 
10""* (blue dash-dotted line) and L{t) = 10"® (black solid line). The two curves are indistinguishable. 
Red dashed line is the rescaled one-dimensional ground state |i?B,iD(2;)|. 



2.2 Simulations 

The radially-symmetric BNLS (|T0)) with d = 2 and = 4 was solved with the initial condition ipo{i") = 2e~^'^~^^ . 
The simulation was run up to L = 0(10^®). Similar results were obtained with d = 3 and cr = 4 (data not 
shown) . 

We next test each item of Conjecture H] numerically: 

1. The position of maximal amplitude rmax(^) — argmax^ {ipl approaches a positive constant as L — >■ 0, 
see Figure [SXl indicating that the solution collapses as a standing ring. 

2. The solution profiles, at the focusing levels of L = 10"* and L — 10~®, rescaled according to 

Acsc.ic4t.p) = L''^"it)iPit,r^,^it)+p-L), L(t) - ll^ll""/', (26) 

are almost indistinguishable, see Figure I5B1 indicating that the collapsing core is self-similar according 
to dill). 

3. Figure [5B] also shows that the self-similar profile of the standing-ring solution is given by i?B.iD(0; the 
one-dimensional ground-state of equation (jl7bl) . 

4. To calculate the blowup rate of ^p, we first assume that L{t) ^ k{Tc — t)P, and find the best fitting k 
and p, see Figure I6AI In this case p « 0.2523, indicating that the blowup rate is a quartic-root or 
slightly faster. 

5. In order to check whether the blowup rate of L is slightly faster than a quartic-root, we compute the 
limit lim L^Lt- Recall that for a quartic-root blowup rate L(t) ~ k{Tc — t)^/* with k > 0, 

lim L^Lt = — - < 0, 

t-i-To 4 

while for a faster-than-a quartic-root blowup rate, see p^ . L^Lf goes to zero. Figure l6Bl shows that 
L^Lt does not approach a negative constant, but increases slowly towards 0~, implying that the blowup 
rate is slightly faster than a quartic root. 
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Figure 6: Blowup rate of the solution of FigureEl A) L as a function of (Tc — <) on a logarithmic scale 
(circles). Solid line is L = 0.774(Tc - t)0-2523 ^ ^^^3 function of 1/L. 



Note that the initial condition ipQ — 2e is quite different from the asymptotic profile tAq^ , indicating 

that the standing-ring ipQ^ profile p5|) is an attractor in the radial case. 



3 Shrinking-ring solutions of the supercritical BNLS (4/<i < cr < 4) 

In this section, we consider the regime 4/(i < cr < 4. In the NLS analogue {2/d < a < 2), the asymptotic 
profile has a "double- lens" radial-phase term, see (jSap . whose explicit form is used in the asymptotic calcu- 
lation of the blowup rate and shrinking rate. In contrast, for the BNLS we do not know the corresponding 
"double-lens" radial phase term, but only the amplitude IV'QbI- Therefore, the results of the asymptotic 
analysis are weaker, and we need to rely on the analogy between the NLS and the BNLS. 



3.1 Informal analysis 

We consider singular ring- type solutions of the supercritical BNLS with A/d < a < A, that undergo a quasi 
self-similar collapse with the asymptotic profile ipQ^^ , whose amplitude is given by 



|V'Qb(^:0 



L^/'^it) 



(27) 



As before, we assume that %p ~ ipg^ in the region Vax — r = 0{L), i.e., for \p\ < Pc ^ We assume 

that a < 1, since otherwise the rings are unstable. Indeed, if a > 1, then p = r/L + o(l), and the rings 
eventually evolve into a peak solution. 
We first derive a lower bound for a: 

Lemma 5. Let A/ d < cr < A, and let ip be a ring-type singular solution of the BNLS equation (|10p . whose 
asymptotic profile is of the form (I27p with a < 1 . Then, 



where 



4-cr 



> 0. 



Therefore, the ring is shrinking, i.e., lim r-^i^^{t) — 0. 
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Proof. First, since 4/c? < ct < 4, then < as < 1- The power of the collapsing core ipq-g is 

f-max + Pc-i(t) 



r - r„ 



L 



2 



'• = I'max-/3c-i(t) 

\Q^{p)\\Lp + roL^f-\Ldp). 

P=-Pa 



In the case a < 1, we have that L\p\ < Lpc <C roL", hence Lp + roL" ^ tqL"'. Therefore, 

II^QbII^ - Li-4/^+"('^"i)(t) . r^^'^'^ r \Q.{p)\'dp. 
In the case a = 1, we have that Lp + roL" = (ro + p)L^ hence 

II^Qb \\l - r |Qe(p)|2 (ro + pf-' dp. 



P=-Pc 



In both cases UqJ\1 = o (ii-4/<T+"(d-i)) . Since UqJ\1 < = Uq\\1 < oo, then L^-i/<^+<^id-i) has 

to be bounded as L — > 0. Therefore, 1 — A/a + a(d — 1) > 0, from which the result follows. □ 



Let 



-Pcoiiapsc = lim inf lim / iV'pr'' ^dr 



be the amount of power that collapses into the singularity. We say that the solution ip undergoes a strong 
collapse if PcoUapsc > 0, and a weak collapse if PcoUapsc = 0. 

Corollary 6. Under the conditions of Lemma \^ 4'Qb undergoes a strong collapse if a — as, and a weak 
collapse if a > Ug. 

Proof. This follows directly from the proof of Lemma [S] □ 



In the NLS with 2/d < a < 2, the shrinking rings undergo a strong collapse with a = a^Ls, see ([Sc]) . 
Therefore, by analogy, we expect that the shrinking rings of the BNLS will also undergo a strong collapse, 
in which case a = a^. 

The blowup rate of singular shrinking-ring solutions of the NLS with 2/d < a <2 in [ FGW07j 

1 1 



P= = 



1 + "nls 2 - 



CT(d-l) 



From the analogy of the BNLS with the NLS (up to the change 2 ^> 4), we expect that the blowup rate of 
singular shrinking rings of the BNLS is 

11 



<j(d-l) 

Therefore, we have the following Conjecture: 



3 + qie 



Conjecture 7. Let d > 1 and A/d < a < A, and let ^ be a singular ring-type solution of the BNLS (jlO[) . 
Then, 

L The solution is quasi self-similar, i.e., "0 ~ ipq^ for r — rmax — 0{L), where ipq^ is given by (j27p . 
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Figure 7: Ring-type singular solution of the supercritical BNLS (fTO|) with d = 2 and a — 8/3. A) rmax 
as a function of the focusing factor 1/L. Solid line is r„iax — 14.9L° '*^^. B) The rescaled solution, 
see p6)). at L{t) = 10~^ (blue dash-dotted line) and L{t) ~ 10^^ (black solid hue). Red dashed line is 
the rescaled one-dimensional ground state |i?B.iD(a^)|- 



2. The solution is a shrinking ring, i.e., lim rniax(0 = 0- 

3. The shrinking rate is 



(28) 



Specifically, < a < 1. 
4-. The blowup rate is 



L{t)^K{n~tY, p^-l-= L_. (29) 



3 + 



Specifically, j < p < ^ . 
In Section 13.21 we provide numerical evidence in support of Conjecture [T] 

3.2 Simulations 

The supercritical BNLS equation with d = 2 and a — 8/3 was solved with the initial condition tpQ — 2e^^^^^^^ . 
The simulation was run up to a focusing level of L{t) — 10^. 
We next test Conjecture [7] numerically, clause by clause. 

1. Figure [7X1 shows that the ring shrinks at a rate of rj^axit) ^ 14.9L"(i) with a sa 0.496, which is close 
to the predicted value of 

4-1 1 



3 



;2-i) 2 



2. In Figure [TB] we plot the solution, rescaled according to (I^Bl) . at the focusing levels 1/L ^ 10 and 1/L = 
100. The two curves are indistinguishable for p = 0(1), but not for all p, showing that the solution 
undergoes a quasi self-similar collapse with the ipq-^ profile (j27p . 
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Figure 8: Blowup rate of the solution of Figure [71 A) Solution width L as a function of (Tc — i) on a 
logarithmic scale (circles). Solid line is L = 0.662(rc — i)"-2844 ^ L'^-^Lt as a function of 1/L. 




Figure 9: The numerical shrinking rate a (circles) for ring-type singular solutions of the BNLS 
with 4/d < cr < 4. The solid line is a = aBicr,d), see A) d = 2 and a = 2, 16/7,8/3, 16/5 and 4. 

B) = 3 and (T = 4/3, 8/5, 2, 8/3 and 4. 



3. Figure [TBI shows that the self-similar profile of the standing-ring solution is close to i?B,iD(0; the 
one-dimcnsional ground-state of equation (jl7b[) . 

4. Figure [83 shows that L{t) - 0.662(rc - t)°-^*'*'*. Therefore, the calculated blowup rate p = 0.282 is 
close to the predicted value of p — g^p^ — 1/3.5 « 0.2857. 



-2.5 1 

c ui .i^ lb ej^cLCLiy pi — ±/o.tJ, cuiiipuLt; iiie iiiiiit liiii ± 

Recall that if L{t) - k{Tc - t)^^^-^, then 



5. In order to check whether the blowup rate of L is exactly p ~ 1/3.5, we compute the limit liin L ' Lt 



lim L^-'^Lt = < 0, 

t-i-Tc 3.5 

while for a faster blowup rate, L^-^Lt — J> 0, and for a slower blowup rate, L^-^Lt — > — oo. Figure [8Bl 
shows that L^'^Lt converges to a negative constant, implying that the blowup rate is exactly p — 1/3.5. 
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0.5 1 0.5 1 

fy fy 

BNLS BNLS 



Figure 10: Same as Figure [HI for the numerical blowup rate p, defined by L{t) ^ c{Tc — t)^, as a 
function of as- Solid line is p = 3^^- The calculated values of p for 1/2 < Ob < 1 are slightly lower 
than the predicted value jrr — ■ 



In Figure |9] we present the numerical values of the shrinking parameter a, defined by r^^axit) ^ cL°' 
for ten different values of (cr, d). In all cases, the value of a is very close to as, see ([28l) . In Figure fTOl 
we present the numerical values of the blowup rate p for the same simulations, and find that they are 
close to p = 1/ (3 + as) for > 0. At = 0+, p has a jump discontinuity to 1/4, in accordance with 
Conjecture [4l The discontinuity 'At — Q [a ~ 4) is a manifestation of the phase transition from shrinking 
to standing rings, see Section FOl 



4 Equal-rate shrinking rings (critical BNLS) 
4.1 Informal Analysis 

We consider singular ring- type solutions of the critical BNLS, that undergo a quasi self-similar collapse with 
the asymptotic profile 

^QBft,r) ^ ;^Qb [ "^l^^^f"^ ) e^'^'\ W^roX". (30) 

Lemma 8. Let ^pQg(t,r), see (j30p . be the asymptotic profile of singular ring-type solutions of the critical 
BNLS (HH). Then, a=l. 

Proof. Theorem [3] implies that the collapsing core of singular solutions of the critical BNLS is self-similar 
in r/L, i.e.. 



m,r)\ ^ 



Therefore, a singular solution of the critical BNLS is ring- type if and only if |^'(p)| attains its maximum at 
some pmax > 0. Hence, rmax(i) = Pmax • Lit). Therefore, a = 1. □ 

By Theorem [21 the blowup rate is lower-bounded by 1/4. We recall that ring- type singular solutions 
of the critical NLS have a square-root blowup rate [F GW05) . Therefore, we expect that ring- type singular 
solutions of the critical BNLS have a quartic-root blowup rate. 

In summary, we conjecture the following: 

Conjecture 9. Let ^ be ring-type singular solution of the critical BNLS (jl2p . Then, 
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-5 5 

p=(r-r )/L 

^ max' 

Figure 11: Ring- type singular solution of the critical BNLS ([T2l) with d ~ 2. A) rmax as a function 
of 1/L. Solid line is fmax = 79.5L^-°'^ . B) The rescaled solution, see at L{t) = 10'^ (blue dashed 
line) and L{t) = 10"^ (black solid hue). 



L The solution undergoes an equal-rate collapse, i.e., r^stx{t) ^ rQL(t). 

2. The solution undergoes a quasi self-similar collapse with the asymptotic profile 



L'^/^it) 

3. The blowup rate is exactly a quartic root, i.e.. 



Qb 



-ds 



L{t) - K{/Tc-t, K>0. 



(31a) 



(31b) 



4.2 Simulations 

The critical BNLS (fT2|) with d ~ 2 was solved with the initial condition -00 — 2.5e'^''^^*'^ . The simulation 
was run up to a focusing level of L{t) = 10^^. We next test Conjecture [9] numerically, clause by clause. 



1. Figure [Tl Al shows that the ring shrinks at a rate of r^nnxit) ^ cL'^{t) with a « 1.02, which is close to 
the predicted value of a = 1. 



2. In Figure HlBI we plot the solution, rescaled according to (l26l) . at the focusing levels 1/L = 10'^ 
and 1/L = 10^. The two curves are indistinguishable, showing that the solution undergoes a quasi 
self-similar collapse with the t/iq^ profile (|31ap . 



3. Figure [T2Al shows that L{t) ~ 0.433(rc — f)^-^^''''. Therefore, the calculated blowup rate is close to a 
quartic root. 

4. By Conjecture m the blowup rate of L{t) should be exactly 1/4, hence L^L{t) should converge to a 
negative constant. However, in Figure [T2B| L^Lt does not converge to a constant, but rather slowly 
decreases away from zero. This indicates that the blowup rate is slower than a quartic-root, which is 
in contradiction with Theorem [2] There are two possible explanations for this: 

(a) It may be that the numerical finding that a is slightly above 1 and p is slightly below 1/4 is an 
artifact of our numerical method, see Section [5j Indeed, in all our simulations for 1/2 < a < 1 in 
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Figure 12: Blowup rate of the solution of Figure [TT] A) Solution width L as a function of (Tc — t) on 
a logarithmic scale (circles). Solid line is L = 0.433(Tc — i)0-2476 g-j j^^^s function of 1/L. 



Figures and [TUl the calculated values of the shrinking rate a were all slightly above as, and the 
blowup rates were slightly below j^^- In those cases, however, these small differences did not 
change the qualitative behavior of the solution. In contrast, a small increment (whether numerical 
or genuine) from a = 1 will drastically change the dynamics, from equal-rate ring-type solutions 
into peak-type solutions. 

(b) It may be that ring-type solutions of the critical BNLS are only meta-stable, having shrinking 
rates a > 1 and a blowup rate slower than 1/4. This does not contradict with Theorem [51 since 
in this case the ring-type solutions will eventually evolve into peak-type solutions with a different 
blowup rate. 

We do not know which of the two options is true. 

5 Numerical Method: Adaptive mesh construction 

In this study, we computed singular solutions of the BNLS equation ([TU]). These solutions become highly- 
localized, so that the spatial scale-difference between the singular region r — rmax ~ 0{L) and the exterior 
regions can be as large as 0{\/L) ^ 10^°. In order to resolve the solution at both the singular and non- 
singular regions, we use an adaptive grid. 

We generate the adaptive grids using the Static Grid Redistribution (SGR) method, which was first intro- 
duced by Ren and Wang [RWOO] , and later simplified and improved by Gavish and Ditkowsky |DG09] . Using 
this approach, the solution is allowed to propagate (self-focus) until it becomes under-resolved. At this stage, 
a new grid, with the same number of grid-points, is generated using De'Boors 'equidistribution principle', 
wherein the grid points {rm} are spaced such that a certain weight function u'i['0] is equidistributed, i.e., 
that 




Wi [ip(r)] dr — const. 



see |RW00llDG09] for details. 
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Algorithm 1 The SGR method, as implemented in [DG09 . 

1. Find a nonlinear coordinate transformation r{x) : [0,1] — )■ [0,R], under which the weight func- 
tion w [ip(r{x))] becomes uniformly distributed. 

2. Transform the solution and equation to the new coordinate system. For example, the second spatial 
derivative in the NLS transforms as 

3. Approximate the equation on a uniform grid {xra\^ using standard finite-differences (or another method 
of choice) . 



The method of |DG09j is given in Algorithm [TJ Note that, since r{x) is nonlinear, the mapping of 
the derivatives of ij) (step 2) involves nonlinear combinations of the derivatives of r. This is not a great 
problem for the NLS, which has only second-order derivatives, but becomes much messier for the biharmonic 
operator (fTTj) . with its many high-order derivatives of ip. 

Therefore, in this study we implement a simplified version of the method of [DG09j . which IS given m 
Algorithm [21 which uses a non-uniform grid in the old-coordinate system, and thereby dispenses with the 
need for transforming the equation, and is much easier to implement in the biharmonic case. 



Algorithm 2 The SGR method, as implemented in this work. 

1. Find a nonlinear coordinate transformation r{x) : [0,1] [0,i?], under which the weight func- 
tion w [^{r{x)y\ becomes uniformly distributed. 

2. Create the uniform grid in the transformed system {xm}- 

3. Create the (highly) non- uniform grid r,„ = f(x„j) in the original (physical) coordinate system. 

4. On the non-uniform grid, approximate the equation using standard (non-uniform) finite-differences. 



We use a third-order accurate finite-difference approximation of the radial biharmonic operator (|lip , with a 
seven-point stencil. 

The method in |DG09j allows control of the fraction of grid points that migrate into the singular region, 
preventing under- resolution at the exterior regions. This is done by using a weight-function W2, which 
penalizes large inter-grid distances. However, we found that this numerical mechanism, while necessary, is 
insufficient for our purposes. In order to understand the reason, let us consider the grid-point spacings Ar„i = 
Tm+i — fm- Using the method of |DG09j with both wi and W2 causes a very sharp bi-partition of the grid 
points - to those inside the singular region, whose spacing is determined by wi and is Ar.,„ — 0{L), and to 
those outside the singular region, whose spacing is determined by W2 and is Ar^ — C(l), see Figure [13AI 
Inside each of these regions, the finite difference approximation we use is well conditioned. However, at the 
transition between these two regions, the finite-difference stencil, seven-points in width, spans grid-spacings 
with 0{1/L) scale-difference — leading to under-resolution which completely violates the validity of the 
finite-difference approximation. 

In order to overcome this limitation, we improve the algorithm of |DG09j by adding a third weight 
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50 100 150 200 250 300 50 100 150 200 250 300 

grid index m grid index m 

(A) old method, L = lO"*^ (B) new method, L = 10" 



Figure 13: The grid-spacing Ar^ obtained using the SGR method of |DG09) for a peak-type singular 
solution of the BNLS. A) The grid generated the original method of |DG09) at focusing level of L = 10^^. 
The Singular and non-singular regions are well-resolved, but the transition region Ar,„ displays a 
discontinuity. At this point, the finite difference operator becomes ill-conditioned. B) same as (A), 
after adding the new penalty function W3, at focusing level L = 10^^^. Even at this much larger 
focusing level, the transition region is now well resolved. 



function 




which penalizes the second-difference A^r„i = Ar,„4.i — Ar,„ operator of the grid locations, allowing for a 
smooth transition between the singular region and the non-singular region, see Fig I13BI 

On the sequence of grids, the equations are solved using a Predictor-Corrector Crank-Nicholson scheme, 
which is second-order in time. 
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